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The energy-level scheme and wave functions of the titanium ions in LaTiOs are calculated using 
crystal-field theory and spin-orbit coupling. The theoretically derived temperature dependence 
and anisotropy of the magnetic susceptibility agree well with experimental data obtained in an 
untwinned single crystal. The refined fitting procedure reveals an almost isotropic molecular field 
and a temperature dependence of the van Vleck susceptibility. The charge distribution of the 3(i- 
electron on the Ti positions and the principle values of the quadrupole moments are derived and 
agree with NMR data and recent measurements of orbital momentum (/) and crystal-field splitting. 
The low value of the ordered moment in the antiferromagnetic phase is discussed. 

PACS numbers: 71.70.Ch, 75.30.Gw, 71.30.+h, 71.27.+a 

Keywords: transit ion-metal oxides, crystal-field analysis, magnetic susceptibility 



I. INTRODUCTION 

In the physics of highly correlated electron systems the 
electronic orbitals and their interactions are in the focus 
of recent experimental and theoretical research, because 
the orbitals play a key role in the coupling of charge and 
spin of the electrons with the lattice. Transition-metal 
oxides, where the shape and anisotropy of the d-electron 
orbitals determine the fundamental electronic properties, 
provide a rich field for this kind of investigations. For ex- 
ample the perovskite titanates ATiOa (with A = Y, La, 
or some trivalent rare-earth ion) are known as realization 
of a Mott insulator. The 3(i^ electronic configuration of 
Ti^+ corresponds to an effectively half-filled conduction 
band, where the on-site Coulomb repulsion inhibits dou- 
ble occupation of the Ti-sites resulting in an insulating 
ground state. ^ Although from their electronic configura- 
tion these titanates seem to be quite simple model sys- 
tems, their orbital properties still have to be resolved 
especially in the case of LaTiOa. 

The debate on the orbital ground state of LaTiOa was 
triggered by its unusual magnetic properties. Below the 
Neel temperature Tn = 146 K,^ LaTiOa reveals a slightly 
canted G-type antiferromagnetic structure with an or- 
dered moment of which is strongly reduced as 
compared to the spin-only value of l/ig and, hence, indi- 
cates a strong importance of the spin-orbit coupling. On 
the other hand the nearly isotropic spin-wave dispersion 
with a small gap of about 3 meV contradicts a dominant 
spin-orbit coupling.^ 

This puzzling situation originates from the fact that 
the orthorhombic GdFeOs structure of LaTiOs deviates 
only weakly from the ideal cubic perovskite structure: 
The quasicubic crystal field of the nearly ideal oxygen 
octahedron surrounding the Ti^+ ion splits the five or- 



bital ?>d levels into a lower t2g triplet and an excited eg 
doublet. The single electron occupies the lower t2g triplet 
and is Jahn- Teller active.^ In principle the Jahn- Teller ef- 
fect is expected to lift the remaining threefold degeneracy 
resulting in a distortion of the oxygen octahedron in favor 
of one of the three orbitals. However, the competing in- 
fluence of spin-orbit coupling cannot be neglected in the 
case of a single electron in a t2g level, as has been outlined 
already by Goodenough'' and by Kugel and Khomskiii^ 
It is important to note that, as long as the orbital triplet 
remains degenerate, the exchange interactions are inher- 
ently frustrated even in a cubic lattice.^ 

To promote possible physics of this degeneracy 
in LaTiOs, an orbital-liquid ground state has been 
suggested. Further detailed theoretical studies^ ^^^^ fa- 
voring the orbital-liquid picture worked out that the frus- 
tration can be resolved via an order-by-disorder mecha- 
nism giving rise to magnetic spin order with disordered 
orbital states. The observed spin- wave excitations were 
found to be in accord with this model. In a differ- 
ent theoretical approach^^'^4ii^*iSii^ the crystal field of 
the La ions caused by the GdFeOs-type distortion has 
been shown to lift the degeneracy of the Ti-t2p-orbitals 
and to stabilize the antiferromagnetic G-type order. In 
the orbital-ground state was derived as ap- 
proximately ^zlii — = {dxy + dyz H- dzx)/ a/3. However, 
Solovye\iii has found that the Hartree-Fock approxima- 
tion alone fails to provide the description of the magnetic 
properties of LaTiOs and YTiOs- 

Several recent experimental investigations strongly 
support the existence of orbital order in LaTiOs- 
Specific-heat, electrical resistivity, thermal-expansion, 
and infrared experiments^^ exhibit anomalies near the 
Neel temperature, which indicate significant structural 
changes and have been interpreted in terms of the 
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influence of orbital order via magneto-elastic interac- 
tions. Transmission-electron microscopy revealed small 
atomic displacements ascribed to a weak Jahn- Teller 
distortion. Detailed x-ray and neutron-diffraction 
studies^^ of crystal and magnetic structure revealed an 
intrinsic distortion of the oxygen octahedra, which leads 
to a large enough splitting of the Ti-t2g triplet state. The 
remeasured magnetic moment ja = 0.57 (5) ja^ turned out 
to be slightly larger than determined before^ The reex- 
amination of the Ti nuclear magnetic resonance spectraiSi 
proves a large nuclear quadrupole splitting, which is as- 
cribed to a rather large quadrupole moment of the 3d 
electrons at the Ti sites. This discarded the earlier 
interpretation^^ of the NMR results in terms of orbital 
degeneracy and clearly favored the orbital order. 

In this communication we perform a detailed anal- 
ysis of the temperature dependence and anisotropy of 
the magnetic susceptibility of LaTiOa, which we ob- 
tained on an untwinned single crystal. In an earlier 
publication^^ it was mentioned that the anisotropy ob- 
served in the paramagnetic regime requires to include 
the spin-orbit coupling into the crystal-field calculation. 
In the present analysis we develop this approach and go 
beyond the Hartree-Fock approximationJI. Besides the 
spin-orbit coupling we are taking into account the Ti-0 
exchange as well. We will show that the obtained orbital- 
order pattern is basically in agreement with NMR data^^ 
and allows to describe consistently the temperature de- 
pendence and anisotropy of the observed experimental 
susceptibility. 



II. CRYSTAL FIELD ANALYSIS 

In LaTiOa the Ti^+ ions (electronic configuration 3d^ , 
spin s = 1/2) are situated in slightly distorted octahedra 
formed by the oxygen ions. The dominant cubic compo- 
nent of the crystal field splits the five Sd-electron states 
into a lower triplet t2g and an upper doublet Cg. The 
low-symmetry component of the crystal field is expected 
to be small with respect to the cubic one and, therefore, 
one may be tempted to analyze the magnetic susceptibil- 
ity using the basis of the t2g states with a fictitious orbital 
momentum / = 1, only^ However, this procedure is not 
convenient for LaTiOa for the following reason: Indeed, 
the wave functions of the fictitious momentum / = 1 are 
defined in a local coordinate system (x, z) with its axes 
parallel to the C4 axes of the non distorted octahedra. 
In the real structure of LaTiOa, there are four different 
fragments TiOe, which are distorted and rotated with 
respect to each other, i.e. the / = 1 basis should be ro- 
tated correspondingly for each of the four inequivalent 
octahedra. During these rotations all Sd-electron states 
are mixed. In this situation it is preferable to stay in the 
crystallographic coordinate system using the full basis of 
3(i-electron states. 

Thus, to determine the energy- level scheme of Ti^+ in 



LaTiOs, we start from the Hamiltonian 

k 

wo = ?(is)+ J2 E^^'^c'W(t?,^). (1) 

The first term denotes the spin-orbit coupling with spin s 
and orbital momentum 1. For Ti^+ the parameter of the 
spin-orbit coupling is expected to be about ^ ^ 200 K.^^i 
The second term represents the crystal field with the 
spherical tensor C^^\i9,(p) = y/27r/{2k ^ l)Yq^^\i^,ip). 
The crystal-field parameters 

= ^a('=)(i?,)(-l)'^ci';)(^„^^.) (2) 

3 

are calculated using available data about the crystal 
structure. ^^-^^-^^ The sum runs over the lattice sites Rj. 

(k) 

The main contributions to the quantities Bq originate 
from the point charges Zj of the lattice and so called 
exchange charges. Hence, the intrinsic parameters of the 
crystal field are given by 

a^'HRj) = -^'^^ai'J{R,). (3) 

The exchange contribution originates from the charge 
transfer from oxygen into the unfilled Sd-shell, i.e. the 
covalence effect, and the direct titanium-oxygen exchange 
coupling:^^-^^ 

= + - ^sh^)- (4) 

Here Ssda , Ssdn 1 and S^s denote the overlap integrals for 
Ti^+(3d^)-02-(2s22/), which are determined in local 
coordinate systems with the z-axis along the titanium- 
oxygen bond. All integrals are calculated using the 
Hartree-Fock wave functions^^ of Ti^+ and 0^~. The 
parameter G = 7.2 is an adjustable parameter, which we 
have extracted from the cubic crystal-field splitting pa- 
rameter lOi^g, which can be assumed as approximately 
similar for all titanium oxides, as e. g. for Ti^+ in AI2O3 
with lODq = 19000 cm-^24 

In LaTiOs there is no inversion symmetry at the oxy- 
gen position and, therefore, each oxygen ion exhibits a 
dipole moment = aE^, where a denotes the polar- 
ization constant^^ and E^ is the electric field of the sur- 
rounding ions at the oxygen site with number i. For 
the oxygen positionsSfiiSi 01 (X = 0.49036, 0.25, Z = 
0.07813) and 02 {x = 0.29144, y = 0.04116, z = 0.71036) 
at T = 298 K, the values of the dipole moments (in 
units of eA) were calculated as dx = —0.093, dy = 0, 
d^ = -0.001 (01) and d^ = 0.036, dy = 0.018, d^ = 0.037 
(02), respectively. The relative signs for the other three 
01 and seven 02 positions change like the signs of the 
corresponding coordinates (X, Z, and x, ^, z), e. g. 
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TABLE I: Contributions to the crystal- field parameters in 
LaTiOs at the Til position |, 0) in units of K 





point charges 


exchange charges 


dipolar 




1527 


720 


-819 




-162-i376 


-301 + i62 


1548 - 2413 


Bi'^ 


-1229 + i2496 


-941 + 2103 


430 + 21525 




-4486 


-7713 


small 


B^^ 


-5828 + ^4105 


10951 + 27733 


small 


Bi'^ 


11325 -^1699 


19452 - 22160 


small 


^(4) 


1827 + i7634 


3407 + 214371 


small 


^(4) 


7638 + il713 


13047 + 22963 


small 



TABLE IL Relative signs of the parameters Bq ^ for Ti2, 
Ti3, and Ti4 with respect to the signs for the Til position 
in LaTiOs 



r(2) 
Bq 


Ti2(0,|,i) 


Ti3(|,0,0) 


Ti4(0,0, 1) 


+ 


+ 


+ 




Re — , Im — 


Re+, Im — 


Re — , Im + 


Bf' 


Re+, Im + 


Re +, Im — 


Re +, Im — 


Bf 


+ 


+ 


+ 




Re — , Im — 


Re+, Im — 


Re — , Im + 


^2 


Re+, Im + 


Re+, Im — 


Re +, Im — 


Bt' 


Re — , Im — 


Re+, Im — 


Re — , Im + 


Bt^ 


Re+, Im + 


Re+, Im — 


Re +, Im — 



for the 01 position (X + 0.5, 0.25, 0.5 - Z) we obtain 
dx = —0.093, dy = 0^ dz = 0.001 and so on. The corre- 
sponding expressions for corrections to the crystal-field 
parameters 5q^^ , B^'' , B^"^ are calculated as usuali^ 

In the crystallographic coordinate system, with the 
Cartesian axes x, and z chosen along the crystal axes 
a = 5.6071 A, 6 = 7.9175 A, and c = 5.6247 A in Pnma 
representation (corresponding to 6, c, and a in Pbnm 
representation, which is used in many papers), respec- 
tively, (values at room temperature 298 K) we obtain 
the crystal-field parameters (in K) for the titanium ion 
in position Til(|, ^,0) as given in Tabled For the other 

(k) 

three titanium positions the absolute values of Bq are 
the same, but their signs are different (cf. Table llT|) . Note 
that the quantum mechanical contributions are compa- 
rable to the classical ones and even dominate for /c = 4. 

Using the crystal-field parameters listed above, for the 
position Til(^, ^,0) we obtain the following energy spec- 
trum of five Kramers doublets with energies £1,2/^6 = 0, 
esA/kB = 2553 K, £5,6/^6 = 3214 K, S7,8/kB = 26773 K, 
and £9,10/^ = 27890 K. This excitation spectrum agrees 
perfectly with results from FIR experiments, which reveal 
a hump in the optical conductivity close to 3000 K.^^ 
It is also in good agreement with the results of recent 
spin-polarized photo-electron spectroscopy experiments, 
which yield a crystal- field splitting of 0.12-0.30 eV, i.e. 
1300-3300 K, of the t2g subshelL^^ The corresponding 
wave functions in \mi,ms) quantization are written as 



TABLE III: Coefficients of the ground- state wave functions 
in LaTiOs at the Til position (|, |, 0) 



i^m^ ,rr 


Is 


ms =T 


ms =i 


mi = 


2 


-0.479 - 20.191 


-0.033 - 20.031 


mi = 


1 


0.136 + 20.025 


0.005 - 2O.O2O 


mi = 





-0.032 + 20.608 


-0.011 +20.030 


mi = 


-1 


0.154 - 20.047 


-0.012 - 20.007 


mi = 


-2 


0.526 - 20.186 


0.048 



follows: 

+2 

1^^) = Yl Yl ^^^ImJ^h^s)- (5) 

mi = -2 nis^hi 

In particular for one of the components of the ground 
doublet of Til (I , ^,0) the coefficients are explicitly given 
in Table ITTTl The other component of the ground state 
can be obtained as Kramers conjugated state. Note that 
the ^-values = 2{ei \ kJz + 2^;, \ Si), = 2 \ {si \ 
kJx + 2sx I £2) I, and gy = 2 \ {ei \ kyly + 2sy \ 82) \ are 
equal for all four titanium positions, i.e. gz — 1.81, gx — 
1.73^ gy = 1.79, where the reduction factors of the orbital 
momentum due to covalency have been assumed as ka = 
1. The relatively small deviation of the g value from the 
spin-only value 2 displays that the orbital momentum is 
rather small, again in agreement with the recent spin- 
resolved photo-emission experiments.^^ 

Figure ^ illustrates the orbital order pattern due to 
the derived ground-state wave function (cf. Table lni|) . 
Basically, this is in agreement with the order patterns 
found by Cwik et al.,^^ by Kiyama and Itoh^ and by 
Pavarini et al}^ However, in those works the wave func- 
tions have been approximated in terms of the t^^g^^ basis 
only, neglecting the spin-orbit coupling. 

Having obtained the orbital ground state, we are able 
to determine the charge distribution at the Ti sites char- 
acterized by the quadrupole moments. The tensor of the 
quadrupole moment per one Ti position is given by 

Qap = ^\e\{r^){3lah-Q5^^). (6) 

Diagonalization of the tensors (5a/3/(|e| (r^)) calculated 
for all four Ti positions yields the same principal values 
equal to Qi = -0.520, Q2 = 0.460, and Qs = 0.060, i. e. 
the charge distribution on the titanium ions is the same in 
the local coordinate systems, which are just rotated with 
respect to each other. The angles of rotations have been 
calculated via the eigenvectors of the tensors Qa(3- The 
components of the unit vectors (n^;, n^, n^) correspond- 
ing to the principal values -0.520 and 0.460 read ni = 
(0.815,0.573,0.086) and n2 = (-0.573,0.746,0.355) at 
the Til site. For -0.520 (0.460) Uy {ux) is reversed at the 
Ti3 and Ti4 sites, whereas is reversed at the Ti2 (Ti2) 
and Ti4 (Ti3) sites. 

It is interesting to know, how the spin is oriented 
with respect to the quadrupole charge distribution. Ac- 
cording to neutron-scattering data^2i24 and susceptibil- 
ity measurements^^ the effective magnetic moment per 
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1^ 




FIG. 1: Orbital order in LaTiOs as derived from the crystal- 
field analysis 

one Ti^+ is about /i^f^ ~ 0.6/i^. The antiferromagneti- 
cally ordered moments are aligned along the c-direction 
and weak ferromagnetism shows up along the 6-direction 
(in Pnma).^^ We suggest that this can be explained 
as follows: Due to the spin-orbit coupling the orienta- 
tions of the titanium magnetic moments are connected 
with the quadrupole ordering. If we assume that the 
spin is aligned perpendicular to the Sd-electron charge- 
distribution plane, i. e. along n2, a ferromagnetic aligne- 
ment along the 6-axis can result from the y component 
of n2, which is positive at all four Ti places, and a G- 
type antiferromagnetic order along the c-axis is favored 
as the sign of the z-component of n2 changes between the 
Ti sites, correspondingly. As neutron scattering detects 
the averaged magnetic moment of the four inequivalent 
Ti places per unit cell with vice versa twisting of the 
quadrupolar moments, the observed fi^^ ~ 0.6/i^ is just 
the projection of the total magnetic moments onto the 
c-direction. 



III. MAGNETIC SUSCEPTIBILITY 

The LaTiOs single crystal, prepared by floating zone 
melting,^ was essentially the same as used previously 
for the thermal-expansion measurements described in 
Ref. The crystallographic axes were determined 

from x-ray Laue pictures. Additional neutron-diffraction 
experiments^^ on the same crystal revealed only a small 
twin-domain of about 5% of the crystal volume, hence 
the crystal can be regarded as practically untwinned. 
The magnetization M{T) was measured in a commercial 



superconducting quantum interference device (SQUID) 
magnetometer (MPMS5, Quantum Design)^ working in 
a temperature range 1.8 < T < 400 K and in magnetic 
fields up to = 50 kOe. 

Figure [21 shows the temperature dependence of the sus- 
ceptibility X = M/H obtained from the LaTiOs single 
crystal in an external field of H = 10 kOe applied along 
the three orthorhombic axes both below Tn (inset) and in 
inverse representation in the paramagnetic regime (main 
frame) . The data have been corrected accounting for the 
diamagnetic background of the sample holder, which was 
measured independently for all three geometries. Below 
the Neel temperature Tn = 146 K, one observes the evo- 
lution of a weak ferromagnetic magnetization of about 
0.02/iB per formula unit with its easy direction along 
the c-axis. The paramagnetic regime is better visible 
in the inverse susceptibility with an approximately linear 
increase above 200 K. Evaluation by a Curie- Weiss be- 
havior A^A/ieff/3^B(T + Ocw), with fil^ = filg^S{S^l) 
yields a Curie- Weiss temperature Ocw ^ 900 K and 
an effective moment /i^^^ ^ 2.6/^3, which is strongly en- 
hanced with respect to the spin-only value of l.TS/ig. For 
an appropriate evaluation we have to take into account 
the preceding energy-level scheme derived from our CF 
analysis. 

Including the external magnetic field, the perturbation 
Hamiltonian is written as 

V -/J^B^aikJa + ^^Sa-faSa) = -HaMa- (7) 

Here the factors fa take into account the molecular field, 
which can be anisotropic for two reasons. The first one 
is because of the anisotropic ^-factors. The second one 
is due to the anisotropy of the effective superexchange 
interaction between the titanium spins, which we take 
in the form X^J^^s^s^. The parameters J^j represent 

a 

the effective superexchange integrals, a = x^y^z. In the 
crystal structure around each Ti^+ ion, there are two 
titanium ions at a distance Ri = 3.958 A, four titanium 
ions at R2 = 3.971 A, and 12 at a distance R3 ~ 5.6 A. 
According to the neutron-scattering data^ ^ ^ 
180 K for ah a. 

The molecular field approximation taking into account 
the six nearest neighbors at distances Ri and R2 yields 

r 6 J^(5q,) (A^q^/q; + 25q,) Cq, , . 

" kBT + QJiSa)^ ~ T + K' ^ ' 

Note that in this approximation the ratios Ca/Qa are 
independent on the exchange coupling Ja and directly 
determined by the spin- and orbital state as 

Cq ^ {ei I kJa + 2Sa I gi) 
Qa {ei I Sa I ^l) 

The ratios Coc/e^ ^ 1-92 and Cy/Sy ^ 1.85, and 
Cz/Qz ^ 1.80, as calculated from the ground state as- 
suming /cq; = 1, indicate again an only small contribu- 
tion of the orbital momentum 1^ to the magnetic suscep- 
tibility. For zero orbital momentum one would obtain 
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T ' 1 ' r 




I I I I I I 

200 300 400 

T(K) 

FIG. 2: Temperature dependence of the inverse susceptibility 
l/x{T) (Inset: x{T) at low temperatures) in LaTiOs for an 
external field oi H — 10 kOe applied along the three crystalo- 
graphic axes a, 6, and c (Pnma). The fits indicated by solid, 
dashed, and dotted lines are described in the text (color on- 
line) . 

For a = z the paramagnetic part of the susceptibility 
can be written as 

Xllra = ^Y.^ei I I si)^ eM-ei/kBT) (10) 
I 

where Z = /cBT^exp(— s^/ZcbT). The van Vleck like 

I 

contribution reads 



The cases a = x^y can be written analogously. In 
addition, one has to take into account the diamag- 
netic susceptibility. It can be estimated from the ionic 
susceptibilities^^ (given in 10~^ emu/mol) of Sr^+ (-15), 
which is isoelectronic to La^+, Ti^+ (-9), and (-12) 
as Xdia — ~^ ' 10~^ emu/mol. 

In Fig. 121 the theoretical description of the data has 
been performed in three steps, as illustrated by the three 
groups of dashed, dotted, and solid lines, respectively. In 
the first step (dashed lines) the exchange coupling is as- 
sumed to be isotropic and used as the only fit parameter 
Joe = J' The reduction factors have been kept fixed at 
/Cq, = 1. With J = 200 K in good agreement with the 
results of neutron scattering, one achieves already a rea- 
sonable description of the susceptibility. It is remarkable 



that absolute value and anisotropy are very well repro- 
duced by this straight-forward calculation. 

In the second step, we allowed a variation of the cova- 
lency parameters ka- With the same exchange constant 
of 200 K and fc^ = 1, /c^ = 0.88, and = 0.95 (dotted 
lines) the description of the relative splitting of the sus- 
ceptibilities between the different axes is improved, but 
the curvature is still not reproduced. Nevertheless, the 
obtained covalency parameters match the values typically 
observed for Ti^+ ions.^^ The resulting ratios Ca/Qa 
change only slightly Cx/Qx ^ 1-92 and Cy/Qy ~ 1.87, 
and Cz/Qz ^ 1-81 with respect to k^ = 1. 

Finally in the third step, the solid lines show the 
fit of the experimental data using in addition the val- 
ues Ca and &a as adjustable parameters. From fit- 
ting we have got Cx/Ox = 2.45, Cy/Qy = 2.29, and 
Cz/Qz = 2.14. These deviations from the nearest- 
neighbor isotropic molecular-field results can be consid- 
ered as a hint for an spin-orbit dependent exchange like 
J{siSj)ll^l'p between the titanium ions. In principle such 
kind of operators are known and have been discussed in 
a number of papers^^'^^'^^'^^ in application to the suscep- 
tibility of the dimer [Ti2Cl9]~^ and a-priory cannot be 
discarded for LaTiOa. Another influence, which to our 
opinion cannot be excluded as well, is the next nearest 
neighbor interaction between the titanium ions. Obvi- 
ously this question should be addressed to further anal- 
ysis, when more experimental information will be ob- 
tained. However, both types of mentioned interactions 
can produce the corrections of a few percent, but we be- 
lieve that the essential physics of the temperature depen- 
dence of magnetic susceptibility and orbital ordering will 
be the same as described above. 

Note that the factor is quite large and, therefore, 
according to Eq. ^2 Xvv dependent on temperature. 
This fact to our knowledge has not been pointed out in 
literature. We think that this situation should be quite 
general for other titanium compounds as well as for vana- 
dium oxides. 



IV. CONCLUSIONS 

In summary the energy splitting and wave functions 
of the Ti^+ Sd^-electron state have been calculated for 
LaTiOs due to the crystal field including spin-orbit cou- 
pling and Ti-0 exchange. From the derived orbital 
ground state we have estimated the quadrupole moments 
at the Ti sites and have deduced the charge-distribution 
picture for the Sd-electrons in the crystallographic coor- 
dinate system. Based on the orientation of the quadrupo- 
lar tensor, it is possible to suggest an explanation for the 
low value of the ordered moment, observed in the antifer- 
romagnetic state. The straight-forward calculation of the 
paramagnetic susceptibility yields the correct anisotropy, 
which we measured in an untwinned LaTiOs single crys- 
tal. 
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